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Abstract 

A general phenomenological theory is presented for the phase behavior of ferromagnetic su- 
perconductors with spin-triplet electron Cooper pairing. The theory describes in details the 
temperature-pressure phase diagrams of real intermetallic compounds exhibiting the remark- 
able phenomenon of coexistence of spontaneous magnetic moment of the itinerant electrons and 
spin-triplet superconductivity. The quantum phase transitions which may occur in these sys- 
tems are also described. The theory allows for a classification of these itinerant ferromagnetic 
superconductors in two types: type I and type II. The classification is based on quantitative 
criteria. The comparison of theory and experiment is performed and outstanding problems are 
discussed. 



1 Introduction 

In the beginning of this century the unconventional superconductivity of spin-triplet type had been 
experimentally discovered in several itinerant ferromagnets. Since then much experimental and the- 
oretical research on the properties of these systems has been accomplished. Here we review the 
phenomenological theory of ferromagnetic unconventional superconductors with spin-triplet Cooper 
pairing of electrons. Some theoretical aspects of the description of the phases and the phase tran- 
sitions in these interesting systems, including the remarkable phenomenon of coexistence of super- 
conductivity and ferromagnetism are discussed with an emphasis on the comparison of theoretical 
results with experimental data. 

The spin-triplet or p-wave pairing allows parallel spin orientation of the fermion Cooper pairs in 
superfluid 3 He and unconventional superconductors pp. For this reason the resulting unconventional 
superconductivity is robust with respect to effects of external magnetic field and spontaneous ferro- 
magnetic ordering, so it may coexist with the latter. This general argument implies that there could 
be metallic compounds and alloys, for which the coexistence of spin-triplet superconductivity and 
ferromagnetism may be observed. 

Particularly, both superconductivity and itinerant ferromagnetic orders can be created by the same 
band electrons in the metal, which means that spin-1 electron Cooper pairs participate in the forma- 
tion of the itinerant ferromagnetic order. Moreover, under certain conditions the superconductivity 
is enhanced rather than depressed by the uniform ferromagnetic order that can generate it, even 
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in cases when the superconductivity does not appear in a pure form as a net result of indirect 
electron-electron coupling. 

The coexistence of superconductivity and ferromagnetism as a result of collective behavior of f- 
band electrons has been found experimentally for some Uranium-based intermetallic compounds as, 
UGe 2 [21 El SI [5], URhGe UCoGe Q5], and Ulr [HI EE]. At low temperature (T ~ 

1 K) all these compounds exhibit thermodynamically stable phase of coexistence of spin-triplet 
superconductivity and itinerant (/-band) electron ferromagnetism (in short, FS phase). In UGe2 
and Ulr the FS phase appears at high pressure (P ~ 1 GPa) whereas in URhGe and UCoGe, the 
coexistence phase persists up to ambient pressure (10 5 Pa = lbar). 

Experiments, carried out in ZrZn 2 [13], also indicated the appearance of FS phase at T < 1 K in 
a wide range of pressures (0 < P ~ 21 kbar). In Zr-based compounds the ferromagnetism and the 
p-wave superconductivity occur as a result of the collective behavior of the <i-band electrons. Later 
experimental results [HI [15] had imposed the conclusion that bulk superconductivity is lacking in 
ZrZn 2 , but the occurrence of a surface FS phase at surfaces with higher Zr content than that in ZrZn 2 
has been reliably demonstrated. Thus the problem for the coexistence of bulk superconductivity with 
ferromagnetism in ZrZn 2 is still unresolved. This raises the question whether the FS phase in ZrZn 2 
should be studied by surface thermodynamics methods or should it be investigated by considering that 
bulk and surface thermodynamic phenomena can be treated on the same footing. Taking into account 
the mentioned experimental results for ZrZn 2 and their interpretation by the experimentalists [T3| [HI 
IT5] we assume that the unified thermodynamic approach can be applied. As an argument supporting 
this point of view let us mention that the spin-triplet superconductivity occurs not only in bulk 
materials but also in quasi-two-dimensional (2D) systems - thin films and surfaces and quasi-lD 
wires (see, e.g., Refs. [IS])- In ZrZn 2 and UGe 2 both ferromagnetic and superconducting orders 
vanish at the same critical pressure P c , a fact implying that the respective order parameter fields 
strongly depend on each other and should be studied on the same thermodynamic basis [T7] . 

Fig. [1] illustrates the shape of the T — P phase diagrams of real intermetallic compounds. The phase 
transition from the normal (N) to the ferromagnetic phase (FM) (in short, N-FM transition) is shown 
by the line Tp(P). The line Tps{P) of the phase transition from FM to FS (FM-FS transition) may 
have two or more distinct shapes. Beginning from the maximal (critical) pressure P c , this line may 
extend, like in ZrZn 2 , to all pressures P < P c , including the ambient pressure P a ; see the almost 
straight line containing the point 3 in Fig. [TJ A second possible form of this line, as known, for 
example, from UGe 2 experiments, is shown in Fig. [TJ by the curve which begins at P ~ P c , passes 
through the point 2, and terminates at some pressure Pi > P a , where the superconductivity vanishes. 
These are two qualitatively different physical pictures: (a) when the superconductivity survives up 
to ambient pressure (type I), and (b) when the superconducting states are possible only at relatively 
high pressure (for UGe 2 , Pi ~ 1 GPa); type II. At the tricritical points 1, 2 and 3 the order of the 
phase transitions changes from second order (solid lines) to first order (dashed lines). It should be 
emphasized that in all compounds, mentioned above, Tpg(P) is much lower than Tp(P) when the 
pressure P is considerably below the critical pressure P c (for experimental data, see Sec. 8). 

In Fig. [TJ the circle C denotes a narrow domain around P c at relatively low temperatures (T < 
300 mK), where the experimental data are quite few and the predictions about the shape of the 
phase transition are not reliable. It could be assumed, as in the most part of the experimental 
papers, that (T = 0, P = P c ) is the zero temperature point at which both lines T F (P) and T F s(P) 
terminate. A second possibility is that these lines may join in a single (N-FS) phase transition 
line at some point (T > 0, P' c < P c ) above the absolute zero. In this second variant, a direct 
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Figure 1: An illustration of T — P phase diagram of p-wave ferromagnetic superconductors (details are 
omitted): N - normal phase, FM - ferromagnetic phase, FS - phase of coexistence of ferromagnetic order 
and superconductivity, Tp{P) and Tps(P) are the respective phase transition lines: solid lines correspond 
to second order phase transitions, dashed lines stand for first order phase transition; 1 and 2 are tricritical 
points; P c is the critical pressure, and the circle C surrounds a relatively small domain of high pressure and 
low temperature, where the phase diagram may have several forms depending on the particular substance. 
The line of the FM-FS phase transition may extend up to ambient pressure (type I ferromagnetic supercon- 
ductors), or, may terminate at T = at some high pressure P = P\ (type II ferromagnetic superconductors, 
as indicated in the figure). 



N-FS phase transition occurs, although this option exists in a very small domain of temperature 
and pressure variations: from point (0, P c ) to point (T > 0,P C ' < P c ). A third variant is related 
with the possible splitting of the point (0,P C ), so that the N-FM line terminates at (0,P C ), whereas 
the FM-FS line terminates at another zero temperature point (0,Pq c ); Pq c < P c . In this case, 
the p-wave ferromagnetic superconductor has three points of quantum (zero temperature) phase 
transitions [HI [19] . 

These and other possible shapes of T — P phase diagrams are described within the framework of 
the general theory of Ginzburg-Landau (GL) type [HI EEl ED] in a conformity with the experimental 
data; see also Ref. |21j. The same theory has been confirmed by a microscopic derivation based 
on a microscopic Hamiltonian including a spin-generalized BCS term and an additional Heisenberg 
exchange term [22]. 

For all compounds, cited above, the FS phase occurs only in the ferromagnetic phase domain of 
the T — P diagram. Particularly at equilibrium, and for given P, the temperature Tp{P) of the 
normal-to-ferromagnetic phase (or N-FM) transition is never lower than the temperature Tps{P) 
of the ferromagnetic-to-FS phase transition (FM-FS transition). This confirms the point of view 
that the superconductivity in these compounds is triggered by the spontaneous magnetization M, in 
analogy with the well-known triggering of the superfluid phase Ax in 3 He at mK temperatures by the 
external magnetic field H. Such "helium analogy" has been used in some theoretical studies (see, 
e.g., Ref. [231 El]), where Ginzburg-Landau (GL) free energy terms, describing the FS phase were 
derived by symmetry group arguments. The non-unitary state, with a non-zero value of the Cooper 
pair magnetic moment, known from the theory of unconventional superconductors and superfluidity 
in 3 He [I], has been suggested firstly in Ref. [23], and later confirmed in other studies [3 El]; recently, 
the same topic was comprehensively discussed in Ref. [25] . 
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For the spin-triplet ferromagnetic superconductors the trigger mechanism was recently examined in 
detail [20], [21] . The system main properties are specified by terms in the GL expansion of form 
Miipjipk, which represent the interaction of the magnetization M = {Mj;j = 1,2,3} with the 
complex superconducting vector field if) = {ipf,j = 1,2,3}. Particularly, these terms are responsible 
for the appearance of superconductivity > 0) for certain T and P values. A similar trigger 
mechanism is familiar in the context of improper ferroelectrics [26] . 

A crucial feature of these systems is the nonzero magnetic moment of the spin-triplet Cooper pairs. 
As mentioned above, the microscopic theory of magnetism and superconductivity in non- Fermi liquids 
of strongly interacting heavy electrons (/ and d band electrons) is either too complex or insufficiently 
developed to describe the complicated behavior in itinerant ferromagnetic compounds. Several au- 
thors (see [2111 Ell 1231 [2U 125]) have explored the phenomeno logical description by a self-consistent 
mean field theory, and here we will essentially use the thermodynamic results, in particular, results 
from the analysis in Refs. [201 EI]- Mean-field microscopic theory of spin-mediated pairing leading 
to the mentioned non-unitary superconductivity state has been developed in Ref. [17] that is in 
conformity with the phenomenological description that we have done. 

The coexistence of s-wave (conventional) superconductivity and ferromagnetic order is a long-standing 
problem in condensed matter physics [23 [2S1 [2S] ■ While the s-state Cooper pairs contain only oppo- 
site electron spins and can easily be destroyed by the spontaneous magnetic moment, the spin-triplet 
Cooper pairs possess quantum states with parallel orientation of the electron spins and therefore can 
survive in the presence of substantial magnetic moments. This is the basic difference in the magnetic 
behavior of conventional (s-state) and spin-triplet superconductivity phases. In contrast to other 
superconducting materials, for example, ternaty and Chevrel phase compounds, where the effect of 
magnetic order on s-wave superconductivity has been intensively studied in the seventies and eighties 
of last century (see, e.g., Refs. (23 123 [21]), in these ferromagnetic compounds the phase transition 
temperature Tp to the ferromagnetic state is much higher than the phase transition temperature 
Tps from ferromagnetic to a (mixed) state of coexistence of ferromagnetism and superconductiv- 
ity. For example, in UGe2 we have Tps ~ 0.8 K versus maximal Tp = 52 K [21 El IH E]. Another 
important difference between the ternary rare earth compounds and the intermetallic compounds 
(UGe2, UCoGe, etc.), which are of interest in this paper, is that the experiments with the latter do 
no give any evidence for the existence of a standard normal-to-superconducting phase transition in 
zero external magnetic field. This is an indication that the (generic) critical temperature T s of the 
pure superconductivity state in these intermetallic compounds is very low (T s <C T^s), if not zero or 
even negative. 

In the reminder of this paper, we present general thermodynamic treatment of systems with itiner- 
ant ferromagnetic order and superconductivity due to spin-triplet Cooper pairing of the same band 
electrons, which are responsible for the spontaneous magnetic moment. The usual Ginzburg-Landau 
(GL) theory of superconductors has been completed to include the complexity of the vector order 
parameter the magnetization M and new relevant energy terms [20, 21J. We outline the T — P 
phase diagrams of ferromagnetic spin-triplet superconductors and demonstrate that in these mate- 
rials two contrasting types of thermodynamic behavior are possible. The present phenomenological 
approach includes both mean- field and spin-fluctuation theory (SFT), as the arguments in Ref. [30] . 
We propose a simple, yet comprehensive, modeling of P dependence of the free energy parameters, 
resulting in a very good compliance of our theoretical predictions for the shape the T — P phase 
diagrams with the experimental data (for some preliminary results, see Ref. [T8"| [TU]). 

The theoretical analysis is done by the standard methods of phase transition theory [31] . Treatment of 
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fluctuation effects and quantum correlations [2U E2] is not included in this study. But the parameters 
of the generalized GL free energy may be considered either in mean-field approximation as here, or as 
phenomenologically renormalized parameters which are affected by additional physical phenomena, 
as for example, spin fluctuations. 

We demonstrate with the help of present theory that we can outline different possible topologies 
for the T — P phase diagram, depending on the values of Landau parameters, derived from the 
existing experimental data. We show that for spin-triplet ferromagnetic superconductors there exist 
two distinct types of behavior, which we denote as Zr-type (or, alternatively, type I) and U-type (or, 
type II); see Fig. [U This classification of the FS, first mentioned in Ref. [18J, is based on the reliable 
interrelationship between a quantitative criterion derived by us and the thermodynamic properties of 
the ferromagnetic spin-triplet superconductors. Our approach can be also applied to URhGe, UCoGe, 
and Ulr. The results shed light on the problems connected with the order of the quantum phase 
transitions at ultra-low and zero temperatures. They also raise the question for further experimental 
investigations of the detailed structure of the phase diagrams in the high-P/low-T region. 



2 Theoretical framework 



Consider the GL free energy functional of the form 



F(if>,B) 



dx 



B 



/sW + /f(M) + A(V, M)+=-- B.M 

07T 



where the fields if), M, and B are supposed to depend on the spatial vector x G V in the volume V 
of the superconductor. In Eq. p]), the free energy density generated by the generic superconducting 
subsystem (if)) is given by 



-M 4 + fi</>T + f ;>>/> 



with 



+K 3 (A^)*(A^), 



(2) 



(3) 



where a summation over the indices is assumed, the symbol Dj = (hd/idxj + 2\e\Aj/c) of 

covariant differentiation is introduced, and Kj are material parameters pp. The free energy density 
/p(JVf) of a standard ferromagnetic phase transition of second order [31] is 



/ F (M) =c f J2 I VMj| 2 + a f M 2 + — M 4 



(4) 



with Cf,bf > 0, and a/ = a(T — Tf), where a>f > and Tf is the critical temperature, corresponding 
of the generic ferromagnetic phase transition. Finally, the energy fi(if), M) produced by the possible 
couplings of if) and M is given by 



Mif), M) = t l0 M.(if) x if)*) + 5 M 2 \if)\ 2 



(5) 
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where the coupling parameter 70 ~ J depends on the ferromagnetic exchange parameter J > 0, [23J 
[21] and 5q is the standard M — 1/) coupling parameter, known from the theory of multicritical 
phenomena [31] and from studies of coexistence of ferromagnetism and superconductivity in ternary 
compounds j27[ [28] . 

As usual, in Eq. ([2]), a s = (T— T s ), where T s is the critical temperature of the generic superconducting 
transition, b s > 0. The parameters u s and v s and Sq may take some negative values, provided the 
overall stability of the system is preserved. The values of the material parameters fj, — (T s , Tf, a s , atf, 
b s , u s , v s , bj, Kj, 7 and S ) depend on the choice of the substance and on intensive thermodynamic 
parameters, such as the temperature T and the pressure P. From a microscopic point of view, the 
parameters /1 depend on the density of states Uf^f) on the Fermi surface. On the other hand 
Up varies with T and P. Thus the relationships (T, P) Up <^ fi, i-e., the functional relations 
/i[Z7jr(T, P)], are of essential interest. While these relations are unknown, one may suppose some 
direct dependence fi(T,P). The latter should correspond to the experimental data. 

The free energy ([1]) is quite general. It has been deduced by several reliable arguments. In order 
to construct Eq. ([I])-© we have used the standard GL theory of superconductors and the phase 
transition theory with an account of the relevant anisotropy of the p-wave Cooper pairs and the 
crystal anisotropy, described by the u s - and t> s -terms in Eq. ([2]), respectively. Besides, we have 
used the general case of cubic anisotropy, when all three components ipj of 1/? are relevant. Note, 
that in certain real cases, for example, in UGe2, the crystal symmetry is tetragonal, t/> effectively 
behaves as a two-component vector and this leads to a considerable simplification of the theory. As 
shown in Ref. [20], the mentioned anisotropy terms are not essential in the description of the main 
thermodynamic properties, including the shape of the T — P phase diagram. For this reason we shall 
often ignore the respective terms in Eq. ([2]). The 70-term triggers the superconductivity (TVf-triger 
effect [20l[2T]) while the 5oM 2 \-ip\ 2 -term makes the model more realistic for large values of M. This 
allows for an extension of the domain of the stable ferromagnetic order up to zero temperatures for 
a wide range of values of the material parameters and the pressure P. Such a picture corresponds 
to the real situation in ferromagnetic compounds ]20] . 

The total free energy ([T]) is difficult for a theoretical investigation. The various vortex and uniform 
phases described by this complex model cannot be investigated within a single calculation but rather 
one should focus on particular problems. In Ref. [23] the vortex phase was discussed with the help of 
the criterion [33] for a stability of this state near the phase transition line T C 2(B), ; see also, Ref. [31]. 
The phase transition line T c2 (H) of a usual superconductor in external magnetic field H = \H\ is 
located above the phase transition line T s of the uniform (Meissner) phase. The reason is that T s is 
defined by the equation a s (T) = 0, whereas T c2 (H) is a solution of the equation \a s \ — HbH, where 
fiB = \e\h/2mc is the Bohr magneton [34]. For ferromagnetic superconductors, where M > 0, one 
should use the magnetic induction B rather than H. In case of H = one should apply the same 
criterion with respect to the magnetization M for small values of near the phase transition line 
T c2 (M); M = \M\. For this reason we should use the diagonal quadratic form [33] corresponding to 
the entire i/> 2 -part of the total free energy functional ([T]). The lowest energy term in this diagonal 
quadratic part contains a coefficient a of the form a = (a s — 70M — 5M 2 ) [35J. Now the equation 
a(T) = defines the critical temperature of the Meissner phase and the equation |a s | = /j^M stands 
for T C 2(M). It is readily seen that these two equations can be written in the same form, provided 
the parameter 7 in a is substituted by 7q = (70 — Thus the phase transition line corresponding 
to the vortex phase, described by the model ([T]) at zero external magnetic field and generated by 
the magnetization M, can be obtained from the phase transition line corresponding to the uniform 
superconducting phase by an effective change of the value of the parameter 70. Both lines have the 
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same shape and this is a particular property of the present model. The variation of the parameter 
70 generates a family of lines. 

Now we propose a possible way of theoretical treatment of the Tps{P) line of the FM-FS phase 
transition, shown in Fig. (1). This is a crucial point in our theory. The phase transition line of the 
uniform superconducting phase can be calculated within the thermodynamic analysis of the uniform 
phases, described by the free energy (JTJ. This analysis is done in a simple variant of the free energy 
(JTJ in which the fields ifi and M do not depend on the spatial vector x. The accomplishment of 
such analysis will give a formula for the phase transition line Tps(P) which corresponds a Meissner 
phase coexisting with the ferromagnetic order. The theoretical result for Tps(P) will contain a 
unspecified parameter 7 . If the theoretical line Tps(P) is fitted to the experimental data for the 
FM-FS transition line corresponding to a particular compound, the two curves will coincide for 
some value of 70, irrespectively on the structure of the FS phase. If the FS phase contains a vortex 
superconductivity the fitting parameter 7o(e//) should be interpreted as 7q but if the FS phase contains 
Meissner superconductivity, 7o( e //) should be identified as 70. These arguments justify our approach 
to the investigation of the experimental data for the phase diagrams of intermetallic compounds with 
FM and FS phases. In the remainder of this paper, we shall investigate uniform phases. 

3 Model considerations 

In the previous section we have justified a thermodynamic analysis of the free energy ([T]) in terms 
of uniform order parameters. Neglecting the ^-dependence of tj> and M, the free energy per unit 
volume, F/V = f(ip, M) in zero external magnetic field (H = 0), can be written in the form 

m,M) = a s ^| 2 + ||^| 4 +y|^T + |El^ + «/M 2 + |M 4 (6) 

+ z 7o M- (t/> x V*) + 5 M 2 M 2 - 

Here we slightly modify the parameter aj by choosing a/ = otf[T n — Tj?(P)], where n — 1 gives 
the standard form of a/, and n = 2 applies for SFT [30] and the Stoner-Wohlfarth model |36j . 
Previous studies [20] have shown that the anisotropy represented by the u s and v s terms in Eq. ([6]) 
slightly perturbs the size and shape of the stability domains of the phases, while similar effects can 
be achieved by varying the b s factor in the b s \t{)\ 4 term. For these reasons, in the present analysis we 
ignore the anisotropy terms, setting u s = v s = 0, and consider b s = b > as an effective parameter. 
Then, without loss of generality, we are free to choose the magnetization vector to have the form 
M = (0,0, M). 

According to the microscopic theory of band magnetism and superconductivity the macroscopic ma- 
terial parameters in Eq. (El) depend in a quite complex way on the density of states at the Fermi 
level and related microscopic quantities [37]. That is why we can hardly use the microscopic charac- 
teristics of these complex metallic compounds in order to elucidate their thermodynamic properties, 
in particular, in outlining their phase diagrams in some details. However, some microscopic sim- 
ple microscopic models reveal useful results, for example, the zero temperature Stoner-type model 
employed in Ref. [35] . 

We redefine for convenience the free energy (jH]) in a dimensionless form by / = //(6/Mq), where 
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M = [otfTJ /bf\ l l 2 > is the value of the magnetization M corresponding to the pure magnetic 
subsystem (-0 = 0) at T = P = and Tf = 7/(0). The order parameters assume the scaling 
m = M/Mq and <p = i/>/[(&//6) 1//4 Mo], and as a result, the free energy becomes 

f = rep 2 + — + tm 2 + — + 27m0i0 2 sin6 l + ^m 2 ^ 2 , (7) 

where 4>j = \<pj\, (ft = \<p\, and 9 = (9 2 — #i) is the phase angle between the complex ip x = 4>ie l61 
and y?2 = <p2^ 2 - Note that the phase angle 0$, corresponding to the third complex field component 
^3 = 4>3^ te3 does not enter explicitly in the free energy /, given by Eq. (J7J), which is a natural result 
of the continuous space degeneration. The dimensionless parameters t, r, 7 and 71 in Eq. (J7J) are 
given by 

t = f n -fj{P), r = K(f-f 9 ), (8) 

where k = a s b 1 ^ 2 /afb 1 ^ 2 TjQ~ 1 , 7 = 7o/[a/Tj 6] 1 / 2 , and 71 = Sq/ '(bbj) 1 ^ 2 . The reduced temperatures 
are f = T/T f0 , f f (P) = T f (P)/T f0 , and f s (P) = T s (P)/T f0 . 

The analysis involves making simple assumptions for the P dependence of the t, r, 7, and 71 pa- 
rameters in Eq. ([7j). Specifically, we assume that only Tf has a significant P dependence, described 
by 



f / (P) = (l-P) 1 /™, (9) 

where P = P/Po and Pq is a characteristic pressure deduced later. In ZrZn 2 and UGe2 the Pq values 
are very close to the critical pressure P c at which both the ferromagnetic and superconducting orders 
vanish, but in other systems this is not necessarily the case. As we will discuss, the nonlinearity 
(n = 2) of Tf(P) in ZrZn 2 and UGe2 is relevant at relatively high P, at which the N-FM transition 
temperature Tp{P) may not coincide with Tf(P); Tp(P) is the actual line of the N-FM phase 
transition, as shown in Fig. (1). The form (Q of the model function Tf(P) is consistent with 
preceding experimental and theoretical investigations of the N-FM phase transition in ZrZn 2 and 
UGe 2 (see, e.g., Refs. [H [2H EH])- Here we consider only non- negative values of the pressure P (for 
effects at P < 0, see, e.g., Ref. [33]). 

The model function (Q is defined for P < Pq, in particular, for the case of n > 1, but we should 
have in mind that, in fact, the thermodynamic analysis of Eq. ((Tj) includes the parameter t rather 
than Tf(P). This parameter is given by 

t(T, P) = f n - 1 + P, (10) 

and is well defined for any P. This allows for the consideration of pressures P > Pq within the free 
energy (j7j). 

The model function Tj(P) can be naturally generalized to Tf(P) = (1 — P^) 1 /" but the present 
needs of interpretation of experimental data do not require such a complex consideration (hereafter 
we use Eq. (J9]) which corresponds to (3 = 1 and a = n). Besides, other analytical forms of P/(P) 
can also be tested in the free energy ([7]), in particular, expansion in powers of P, or, alternatively, 
in (1 — P) which satisfy the conditions T/(0) = 1 and 7/(1) = 0. Note, that in URhGe the slope 
of Tp(P) ~ Tf(P) is positive from P = up to high pressures [5] and for this compound the form 
(JHJ) of Tf(P) is inconvenient. Here we apply the simplest variants of P-dependence, namely, Eqs. ([HD 
and ([TO]). 
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In more general terms, all material parameters (r, t, 7, . . . ) may depend on the pressure. We suppose 
that a suitable choice of the dependence of t on P is enough for describing the main thermodynamic 
properties and this supposition is supported by the final results, presented in the remainder of 
this paper. But in some particular investigations one may need to introduce a suitable pressure 
dependence of other parameters. 



4 Stable phases 



The simplified model ([7]) is capable of describing the main thermodynamic properties of spin-triplet 
ferromagnetic superconductors. For r > 0, i.e., T > T s , there are three stable phases |20j: (i) 
the normal (N-) phase, given by <fi = m = (stability conditions: t > 0, r > 0); (ii) the pure 
ferromagnetic phase (FM phase), given by m = (— t) 1 / 2 > 0, = 0, which exists for t < and is 
stable provided r > and r > (7^ + 7|t| 1//2 ), and (iii) the already mentioned phase of coexistence of 
ferromagnetic order and superconductivity (FS phase), given by sm8 = =pl, 3 = 0, <fii = 4>2 = 4>/V2, 
where 



2 = k(T s — T) ± 7m — 7!m 2 > 0. 
The magnetization m satisfies the equation 



C3777, ± C2?n + C\VTL ± Cq = 



with coefficients Co = 7k(T — f s ), 

ci = 2 



f n + K ll (f s -f)+P-l- 

c 2 = 377i, c 3 = 2(l-7 2 ). 
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(12) 

(13) 
(14) 



Table 1. Theoretical results for the location [(T, P) - reduced coordinates] of the tricritical points A 
= (Ta, Pa) an d B = (Tb,Pb), the critical-end point C = (Tc,Pc), and the point of temperature max- 
imum, max =(T m , P m ) on the curve Tps(P) of the FM-FS phase transitions of first and second orders (for 



details, see Sec. 5). The first column shows T/v = T, 



(A,B,C,m)- 



The second column stands for tjv = t 



(A,B,C,m)- 



The reduced pressure values P(A,B,c,m) °f points A, B, C, and max are denoted by P/v(n): n = 1 stands for 
the linear dependence Tf(P), and n = 2 stands for the nonlinear Tf(P) and t(T), corresponding to SFT. 



N 


T N 


tN 


Pn 


(n) 


A 


T s 




7 2 /2 


1 - 


T™ + 7 2 /2 


B 


T s + 7 2 (2 + 7i )/4k(1- 


H7i) 2 


- 7 2 /4(l + 7i) 2 


1 - 


T™-7 2 /4(l+7i) 2 


C 


f s + 7 2 /4k(1 + 71) 







1 - 




max 


f s + 7 2 /4k 7 i 




- 7 2 /4 7l 2 


1 - 


tL-i 2 /H 



The FS phase contains two thermodynamically equivalent phase domains that can be distinguished 
by the upper and lower signs (±) of some terms in Eqs. (TTTT) and (fl~2l) . The upper sign describes the 
domain (labelled bellow again by FS), where m > 0, sin# = —1, whereas the lower sign describes 
the conjunct domain FS*, where m < and sin# = 1 (for details, see, Ref. [20J). Here we consider 
one of the two thermodynamically equivalent phase domains, namely, the domain FS, which is stable 
for m > (FS* is stable for m < 0). This "one-domain approximation" correctly presents the main 
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thermodynamic properties described by the model ([H]), in particular, in the case of a lack of external 
symmetry breaking fields. The stability conditions for the FS phase domain given by Eqs. fllll) and 
(II2D are 7M > 0, 

k(T s - f) ± 7m - 271m 2 > 0, (15) 

and 

3(1 - 7l 2 )m 2 + 3 77l m + T n - 1 + P + « 7l (T s - T) - j- > 0. (16) 

These results are valid whenever Tf(P) > T S (P), which excludes any pure superconducting phase 
{ip ^ 0, m = 0) in accord with the available experimental data. 

For r < 0, and t > the models (jSJ) and ([7]) exhibit a stable pure superconducting phase (0i = </> 2 = 
m = 0, 0| = — r) [20]. This phase may occur in the temperature domain T/(P) < T < T s . For 
systems, where 7/(0) ^ T s , this is a domain of pressure in a very close vicinity of Po ~ P c , where 
Tp{P) ~ Tf(P) decreases up to values lower than T s . Of course, such a situation is described by 
the model ([7]) only if T s > 0. This case is interesting from the experimental point of view only when 
T s > is enough above zero to enter in the scope of experimentally measurable temperatures. Up 
to date a pure superconducting phase has not been observed within the accuracy of experiments on 
the mentioned metallic compounds. For this reason, in the reminder of this paper we shall often 
assume that the critical temperature T s of the generic superconducting phase transition is either 
non-positive (T s < 0), or, has a small positive value which can be neglected in the analysis of the 
available experimental data. 

The negative values of the critical temperature T s of the generic superconducting phase transition are 
generally possible and produce a variety of phase diagram topologies (Sec. 5). Note, that the value 
of T s depends on the strength of the interaction mediating the formation of the spin-triplet Cooper 
pairs of electrons. Therefore, for the sensitiveness of such electron couplings to the crystal lattice 
properties, the generic critical temperature T s depends on the pressure. This is an effect which might 
be included in our theoretical scheme by introducing some convenient temperature dependence of T s . 
To do this we need information either from experimental data or from a comprehensive microscopic 
theory. 

Usually, T s < is interpreted as a lack of any superconductivity but here the same non-positive 
values of T s are effectively enhanced to positive values by the interaction parameter 7 which triggers 
the superconductivity up to superconducting phase-transition temperatures T FS (P) > 0. This is 
readily seen from Table 1, where we present the reduced critical temperatures on the FM-FS phase 
transition line Tps(P), calculated from the present theory, namely, T m - the maximum of the curve 
Tps(P) (if available, see Sec. 5), the temperatures T4 and Tb, corresponding to the tricritical points 
A = {Ta,Pa) and B = (Tb,Pb), and the temperature To, corresponding to the critical-end point 
C = (T c ,Pc)- The theoretical derivation of the dependence of the multicritical temperatures Ta, 
T b and T c on 7, ji, k, and T s , as well as the dependence of T m on the same model parameters 
is outlined in Sec. 5. All these temperatures as well as the whole phase transition line Tps{P) are 
considerably boosted above T s owing to positive terms of order 7 2 . If T s < 0, the superconductivity 
appears, provided T m > 0, i.e., when 7 2 /4k7i > \T S \ (see Table 1). 
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5 Temperature-pressure phase diagram 



Although the structure of the FS phase is quite complicated, some of the results can be obtained in 
analytical form. A more detailed outline of the phase domains, for example, in T — P phase diagram, 
can be done by using suitable values of the material parameters in the free energy ((7j): P , Ty , T s , 
k, 7, and 71. Here we present some of the analytical results for the phase transition lines and the 
multi-critical points. Typical shapes of phase diagrams derived directly from Eq. ([7]) are given in 
Figs. 2-7. Figure 2 shows the phase diagram calculated from Eq. ([7]) for parameters, corresponding 
to the experimental data [13] for ZrZn2. Figures 3 and 4 show the low-temperature and the high- 
pressure parts of the same phase diagram (see Sec. 7 for details). Figures 5-7 show the phase diagram 
calculated for the experimental data [2], S] of UGe2 (see Sec. 8). In ZrZn 2 , UGe2, as well as in UCoGe 
and Ulr, critical pressure P c exists, where both superconductivity and ferromagnetic orders vanish. 

As in experiments, we find out from our calculation that in the vicinity of Po ~ P c the FM-FS phase 
transition is of fist order, denoted by the solid line BC in Figs. 3, 4, 6, and 7. At lower pressure the 
same phase transition is of second orderq shown by the dotted lines in the same figures. The second 
order phase transition line T F s(P) separating the FM and FS phases is given by the solution of the 
equation 



f FS (P) = f s + ixtpsiP) + ^[-t FS (P)] 1/2 , (17) 

where t FS (P) = t(T FS ,P) < 0, 7 = j/k, 71 = 71/ft, and < P < Pg; Pb is the pressure 
corresponding to the multi-critical point B, where the line T F s(P) terminates, as clearly shown in 
Figs. 4 and 7). Note, that Eq. (|T7j) strictly coincides with the stability condition for the FM phase 
with respect to appearance of FS phase [20] ■ 

Additional information for the shape of this phase transition line can be obtained by the derivative 
p = df FS (P)/dP, namely, 

~_ Ps + ii-l/2(-t F s) 1/2 m 
l-nf n FS l [li-l/2[{-t FS Yl^ 

where p s = df s (P)/dP. Note, that Eq. (PJ is obtained from Eqs. ffTU]) and (JTTj). 

The shape of the line T F s(P) can vary depending on the theory parameters (see, e.g., Figs. 3 and 6). 
For certain ratios of 7, 71, and values of p s , the curve T F s(P) exhibits a maximum T m = T F s(P m ), 
given by p(p s , T m , P m ) = 0. This maximum is clearly seen in Figs. 6 and 7. To locate the maximum 
we need to know p s . We have already assumed T s does not depend on P, as explained above, 
which from the physical point of view means that the function T S (P) is flat enough to allow the 
approximation T s « without a substantial error in the results. From our choice of P-dependence 
of the free energy [Eq. (J7J] parameters, it follow that p s = 0. 

Setting p s = p = in Eq. ( TT8|) we obtain 

7 2 

t(T m , P m ) 7^9; (19) 
471 

namely, the value t m (T, P) = t(T m , P m ) at the maximum T m (P m ) of the curve T F g(P). Substituting 
t m back in Eq. (|17|) we obtain T m , and with its help we also obtain the pressure P m , both given in 
Table 1, respectively. 
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Figure 2: T — P diagram of ZrZn 2 calculated for T s = 0, Tj = 28.5 K, P = 21 kbar 



K = 10, 

7 = 2t~! w 0.2, and n — 1. The dotted line represents the FM-FS transition and the dashed line 
stands for the second order N-FM transition. The dotted line has a zero slope at P = 0. The 
low-temperature and high-pressure domains of the FS phase are seen more clearly in the following 
Figs. 3 and 4. 



We want to draw the attention to a particular feature of the present theory that the coordinates 
T m and P m of the maximum (point max) at the curve Tps(P) as well as the results from various 
calculations with the help of Eqs. (fTT|) and fflBl are expressed in terms of the reduced interaction 
parameters 7 and 71. Thus, using certain experimental data for T m , P m , as well as Eqs. ffTTj) and ([TBI 
for T FS , T s , and the derivative p at particular values of the pressure P, 7 and 71 can be calculated 
without any additional information, for example, for the parameter k. This property of the model 
(J7J) is quite useful in the practical work with the experimental data. 

The conditions for existence of a maximum on the curve T F s{P) can be determined by requiring 
P m > 0, and T m > and using the respective formulae for these quantities, shown in Table 1. This 
max always occurs in systems where Tps(O) < and the low-pressure part of the curve Tps(P) 
terminates at T = for some non- negative critical pressure Pq c (see Sec. 6). But the max may 
occur also for some sets of material parameters, when T FS (0) > (see Fig. 3, where P m = 0). All 
these shapes of the line T FS (P) are described by the model ([7]). Irrespectively of the particular 
shape, the curve T FS (P) given by Eq. (JTTj) always terminates at the tricritical point (labeled B), 
with coordinates (Pb,T b ) (see, e.g., Figs. 4 and 7). 

At pressure P > Pb the FM-FS phase transition is of first order up to the critical-end point C. For 
Pb < P < Pc the FM-FS phase transition is given by the straight line BC (see, e.g., Figs. 4 and 
7). The lines of all three phase transitions, N-FM, N-FS, and FM-FS, terminate at point C. For 
P > Pc the FM-FS phase transition occurs on a rather flat smooth line of equilibrium transition of 
first order up to a second tricritical point A with P4 ~ Pq and Ta ~ 0. Finally, the third transition 
line terminating at the point C describes the second order phase transition N-FM. The reduced 
temperatures Tjv and pressures Pn, N = (A, B, C, max) at the three multi-critical points (A, B, and 
C), and the maximum T m (P m ) are given in Table 1. Note that, for any set of material parameters, 
T A <T C <T B < T m and P m < P B < P c < Pa- 

There are other types of phase diagrams, resulting from model ([7]). For negative values of the generic 
superconducting temperature T s , several other topologies of the T — P diagram can be outlined. The 
results for the multicritical points, presented in Table 1, shows that, when T s lowers below T = 0, 
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Figure 3: Details of Fig. 2 with expanded temperature scale. The points A, B, C are located in the 
high-pressure part (P ~ P c ~ 21 kbar). The max point is at P w kbar. The FS phase domain 
is shaded. The dotted line shows the second order FM-FS phase transition with P m w 0. The solid 
straight line BC shows the fist-order FM-FS transition for P > Pg. The quite flat solid line AC 
shows the first order N-FS transition (the lines BC and AC are more clearly seen in Fig. 4. The 
dashed line stands for the second order N-FM transition. 



T c also decreases, first to zero, and then to negative values. When T c = the direct N-FS phase 
transition of first order disappears and point C becomes a very special zero-temperature multicritical 
point. As seen from Table 1, this happens for T s = — 7 2 T/(0)/4k(1 + 71). The further decrease of T s 
causes point C to fall below the zero temperature and then the zero-temperature phase transition of 
first order near P c splits into two zero-temperature phase transitions: a second order N-FM transition 
and a first order FM-FS transition, provided Tg still remains positive. 

At lower T s also point B falls below T = and the FM-FS phase transition becomes entirely of 
second order. For very extreme negative values of T s , a very large pressure interval below P c may 
occur where the FM phase is stable up to T = 0. Then the line Tps{P) will exist only for relatively 
small pressure values (P <C P c ). This shape of the stability domain of the FS phase is also possible 
in real systems. 



6 Quantum phase transitions 

We have shown that the free energy ([6]) describes zero temperature phase transitions. Usually, 
the properties of these phase transitions essentially depend on the quantum fluctuations of the 
order parameters. For this reason the phase transitions at ultralow and zero temperature are called 
quantum phase transitions [311 [32] . The time-dependent quantum fluctuations (correlations) which 
describe the intrinsic quantum dynamics of spin-triplet ferromagnetic superconductors at ultralow 
temperatures are not included in our consideration but some basic properties of the quantum phase 
transitions can be outlines within the classical limit described by the free energy models and ([7]). 
Let we briefly clarify this point. 

The classical fluctuations are entirely included in the general GL functional ([I])-© but the quantum 
fluctuations should be added in a further generalization of the theory. Generally, both classical 
(thermal) and quantum fluctuations are investigated by the method of the renormalization group 
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Figure 4: High-pressure part of the phase diagram of ZrZn 2 , shown in Fig. 1. The thick solid lines 
AC and BC show the first-order transitions N-FS, and FM-FS, respectively. Other notations are 
explained in Figs. 2 and 3. 



(RG) [31], which is specially intended to treat the generalized action of system, where the order 
parameter fields (cp and M) fluctuate in time t and space x [HI [32]. These effects, which are 
beyond the scope of the paper, lead either to a precise treatment of the narrow critical region 
in a very close vicinity of second order phase transition lines or to a fluctuation-driven change in 
the phase-transition order. But the thermal fluctuations and quantum correlation effects on the 
thermodynamics of a given system can be unambiguously estimated only after the results from 
counterpart simpler theory, where these phenomena are not present, are known and, hence, the 
distinction in the thermodynamic properties predicted by the respective variants of the theory can 
be established. Here we show that the basic low-temperature and ultralow-temperature properties of 
the spin-triplet ferromagnetic superconductors, as given by the preceding experiments, are derived 
from the model without any account of fluctuation phenomena and quantum correlations. The 
latter might be of use in a more detailed consideration of the close vicinity of quantum critical points 
in the phase diagrams of ferromagnetic spin-triplet superconductors. Here we show that the theory 
predicts quantum critical phenomena only for quite particular physical conditions whereas the low- 
temperature and zero-temperature phase transitions of first order are favored by both symmetry 
arguments and detailed thermodynamic analysis. 

There is a number of experimental [9j HQ] and theoretical [TTJ HTJ H2] investigations of the problem 
for quantum phase transitions in unconventional ferromagnetic superconductors, including the men- 
tioned intermetallic compounds. Some of them are based on different theoretical schemes and do not 
refer to the model ([6]). Others, for example, those in Ref. [JT] reported results about the thermal 
and quantum fluctuations described by the model ([6]) before the comprehensive knowledge for the 
results from the basic treatment reported in the present investigation. In such cases one could not 
be sure about the correct interpretation of the results from the RG and the possibilities for their 
application to particular zero-temperature phase transitions. Here we present basic results for the 
zero-temperature phase transitions described by the model 

The RG investigation [?T] has demonstrated up to two loop order of the theory that the thermal 
fluctuations of the order parameter fields rescale the model ([6]) in a way which corresponds to first 
order phase transitions in magnetically anisotropic systems. This result is important for the metallic 
compounds we consider here because in all of them magnetic anisotropy is present. The uniaxial 
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Figure 5: T — P diagram of UGe 2 calculated taking T s = 0, Tf = 52 K, P = 1.6 GPa, k — 4, 
7 = 0.0984, 7i = 0.1678, and n = 1. The dotted line represents the FM-FS transition and the 
dashed line stands for the N-FM transition. The low-temperature and high-pressure domains of the 
FS phase are seen more clearly in the following Figs. 6 and 7. 

magnetic anisotropy in ZrZn 2 is much weaker than in UGe2 but cannot be neglected when fluctuation 
effects are accounted for. Owing to the particular symmetry of model ([6]), for the case of magnetic 
isotropy (Heisenberg symmetry), the RG study reveals an entirely different class of (classical) critical 
behavior. Besides, the different spatial dimensions of the superconducting and magnetic quantum 
fluctuations imply a lack of stable quantum critical behavior even when the system is completely mag- 
netically isotropic. The pointed arguments and preceding results lead to the reliable conclusion that 
the phase transitions, which have already been proven to be first order in the lowest-order approxi- 
mation, where thermal and quantum fluctuations are neglected, will not undergo a fluctuation-driven 
change in the phase transition order from first to second. Such picture is described below, in Sec. 8, 
and it corresponds to the behavior of real compounds. 

Our results definitely show that the quantum phase transition near P c is of first order. This is valid 
for the whole N-FS phase transition below the critical-end point C, as well as the straight line BC. 
The simultaneous effect of thermal and quantum fluctuations do not change the order of the N-FS 
transition, and it is quite unlikely to suppose that thermal fluctuations of the superconductivity field 
if} can ensure a fluctuation-driven change in the order of the FM-FS transition along the line BC. 
Usually, the fluctuations of if} in low temperature superconductors are small and slightly influence 
the phase transition in a very narrow critical region in the vicinity of the phase-transition point. 
This effect is very weak and can hardly be observed in any experiment on low-temperature supercon- 
ductors. Besides, the fluctuations of the magnetic induction B always tend to a fluctuation- induced 
first-order phase transition rather than to the opposite effect - the generation of magnetic fluctuations 
with infinite correlation length at the equilibrium phase-transition point and, hence, a second order 
phase transition J3TJ [43] . Thus we can quire reliably conclude that the first-order phase transitions 
at low-temperatures, represented by the lines BC and AC in vicinity of P c do not change their order 
as a result of thermal and quantum fluctuation fluctuations. 

Quantum critical behavior for continuous phase transitions in spin-triplet ferromagnetic supercon- 
ductors with magnetic anisotropy can therefore be observed at other zero-temperature transitions, 
which may occur in these systems far from the critical pressure P c . This is possible when Tps(0) = 
and the Tps{P) curve terminates at T = at one or two quantum (zero-temperature) critical points: 
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Figure 6: Low-temperature part of the T — P phase diagram of UGe2, shown in Fig. 5. The points 
A, B, C are located in the high-pressure part (P ~ P c ~ 1.6 GPa). The FS phase domain is shaded. 
The thick solid lines AC and BC show the first-order transitions N-FS, and FM-FS, respectively. 
Other notations are explained in Figs. 2 and 3. 



Poc < P m - "lower critical pressure", and Pq c > P m - "upper critical pressure." In order to obtain 
these critical pressures one should solve Eq. (TTTI) with respect to P, provided Tps(P) = 0, T m > and 
P m > 0, namely, when the continuous function Tps{P) exhibits a maximum. The critical pressure 
Pq c is bounded in the relatively narrow interval (P m , Pb) and can appear for some special sets of 
material parameters (r, t, 7, 71). In particular, as our calculations show, Pq c do not exists for T s > 0. 



7 Criteria for type I and type II spin-triplet 
ferromagnetic superconductors 

The analytical calculation of the critical pressures Pq c and Pq c for the general case of T s 7^ leads 
to quite complex conditions for appearance of the second critical field Pq c . The correct treatment 
of the case T s 7^ can be performed within the entire two-domain picture for the phase FS (see, 
also, Ref. [20] )■ The complete study of this case is beyond our aims but here we will illustrate our 
arguments by investigation of the conditions, under which the critical pressure P oc occurs in systems 
with T s « 0. Moreover, we will present the general result for P 0c > and Pq c > in systems where 

Setting Tps(P 0c ) = in Eq. (JTTj) we obtain the following quadratic equation, 

7im^ c - 7m 0c - f s = 0, (20) 

for the reduced magnetization, 

m 0c =[-t(0,P oc )] 1 /2 = (l_i3 0c) i/2 (21) 

and, hence, for P 0c . For T s 7^ 0, Eqs. (I2"U|) and (|2TT) have two solutions with respect to P 0c . For 
some sets of material parameters these solutions satisfy the physical requirements for P 0c and Pq c 
and can be identified with the critical pressures. The conditions for existence of Pq c and Pq c can 
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Figure 7: High-pressure part of the phase diagram of UGe2, shown in Fig. 4. Notations are explained 
in Figs. 2, 3, 5, and 6. 

be obtained either by analytical calculations or by numerical analysis for particular values of the 
material parameters. 

For T s — 0, the trivial solution P 0c = 1 corresponds to Pq c — Pq > Pb and, hence, does not satisfy 
the physical requirements. The second solution, 



0c 



~2 

it 



is positive for 



£>1 

7 



(22) 



(23) 



and, as shown below, it gives the location of the quantum critical point (T = 0,P 0c < P m ). At this 
quantum critical point, the equilibrium magnetization m 0c is given by m 0c = 7/71 and is twice bigger 
that the magnetization m m = 7/271 ([20]) at the maximum of the curve Tps{P)- 

To complete the analysis we must show that the solution (|22|) satisfies the condition Pq c < P m . By 
taking P m from Table 1, we can show that solution (|22|) satisfies the condition P 0c < P m for n — 1, if 



and for n = 2 (SFT case), when 



7i < 3k > 



7 < 2y/3K. 



(24) 



(25) 



Finally, we determine the conditions under which the maximum T m of the curve Tps{P) occurs at 
non- negative pressures. For n = 1, we obtain that P m > for n = 1, if 



7i > 1 f 1+ 7i V /2 
7 ~ 2 V k 



(26) 



whereas for n = 2, the condition is 
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Obviously, the conditions ( I2"3"|) - (j2"7|) are compatible with one another. The condition f l2T)|) is weaker 
than the condition Eq. (|23|) . provided the inequality (|24p is satisfied. The same is valid for the 
condition (|27|) if the inequality ( 125|) is valid. In Sec. 8 we will show that these theoretical predictions 
are confirmed by the experimental data. 

Doing in the same way the analysis of Eq. (fTTj) . some results may easily obtained for T s ^ 0. In this 
more general case the Eq. ( |T7j) has two nontrivial solutions, which yield two possible values of the 
critical pressure 



The relation Pq c (-) > Pqc(+) is always true. Therefore, to have both Poc(±) > 0, it is enough to require 
Poc(+) > 0. Having in mind that for the phase diagram shape, we study T m > 0, and according to the 
result for T m in Table 1, this leads to the inequality T s > — 7 2 /4k7 1 . So, we obtain that Pq c {+) > 
will exist, if 



which generalizes the condition (|23|) . 

Now we can identify the pressure Pq c {+) with the lower critical pressure Po c , and Poc(-) with the 
upper critical pressure Pq c . Therefore, for wide variations in the parameters, theory ([6]) describes a 
quantum critical point P oc , that exists, provided the condition f[2"9~j) is satisfied. The quantum critical 
point (T = 0, Pq c ) exists in UGe2 and, perhaps, in other p-wave ferromagnetic superconductors, for 
example, in Ulr. 

Our results predict the appearance of second critical pressure - the upper critical pressure P' oc that 
exists under more restricted conditions and, hence, can be observed in more particular systems, where 
T s < 0. As mentioned in Sec. 5, for very extreme negative values of T s , when Tb < 0, the upper 
critical pressure Pq c > occurs, whereas the lower critical pressure P$ c > does not appear. Bue 
especially this situation should be investigated in a different way, namely, one should keep Tps(0) 
different from zero in Eq. (fTTj) . and consider a form of the FS phase domain in which the curve 
Tps{P) terminates at T = for Pq c > 0, irrespective of whether the maximum T m exists or not. In 
such geometry of the FS phase domain, the maximum T(P m ) may exist only in quite unusual cases, 
if it exists at all. 

Using criteria like ( 1231) in Sec. 8.4 we classify these superconductors in two types: (i) type I, when the 
condition (|23j) is satisfied, and (ii) type II, when the same condition does not hold. As we show in 
Sec. 8.2, 8.3 and 8.4, the condition (|23l) is satisfied by UGe2 but the same condition fails for ZrZn 2 . 
For this reason the phase diagrams of UGe2 and ZrZn 2 exhibit qualitatively different shapes of the 
curves TpS(P). For UGe2 the line TpS(P) has a maximum at some pressure P > 0, whereas the 
line TpS(P), corresponding to ZrZn2, does not exhibit such maximum (see also Sec. 8). 




(28) 




kT s 



(29) 
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The quantum and thermal fluctuation phenomena in the vicinities of the two critical pressures Pq c 
and Pq c need a nonstandard RG treatment because they are related with the fluctuation behavior of 
the superconducting field ip far below the ferromagnetic phase transitions, where the magnetization 
M does not undergo significant fluctuations and can be considered uniform. The presence of uniform 
magnetization produces couplings of M and if) which are not present in previous RG studies and 
need a special analysis. 

8 Application to metallic compounds 

8.1 Theoretical outline of the phase diagram 

In order to apply the above displayed theoretical calculations, following from free energy ([7]), for the 
outline of T — P diagram of any material, we need information about the values of Pq, Tfo, T s , k 7, 
and 71. The temperature Tfo can be obtained directly from the experimental phase diagrams. The 
pressure Pq is either identical or very close to the critical pressure P c , for which the N-FM phase 
transition line terminates at T ~ 0. The temperature T s of the generic superconducting transition 
is not available from the experiments because, as mentioned above, pure superconducting phase not 
coexisting with ferromagnetism has not been observed. This can be considered as an indication 
that T s is very small and does not produce a measurable effect. So the generic superconducting 
temperature will be estimated on the basis of the following arguments. For Tf(P) > T s we must 
have T S (P) = at P > P c , where Tf(P) < 0, and for < P < Pq, T s < Tc- Therefore for materials 
where Tc is too small to be observed experimentally, T s can be ignored. 

As far as the shape of FM-FS transition line is well described by Eq. (fl7|) . we will make use of 
additional data from available experimental phase diagrams for ferroelectric superconductors. For 
example, in ZrZn 2 these are the observed values of Tps(0) and the slope po = [dTFs(P)/dP]o = 
(Tfo/Po)po at P = 0; see Eq. lITTj) . For UGe2, where a maximum (T m ) is observed on the phase- 
transition line, we can use the experimental values of T m , P m , and Pq c . The interaction parameters 
7 and 71 are derived using Eq. ffTTj) . and the expressions for T m , P m , and p , see Table 1. The 
parameter k is chosen by fitting the expression for the critical-end point Tc- 

8.2 ZrZn 2 

Experiments for ZrZn2 [13J gives the following values: Tfo = 28.5 K, Tps(0) = 0.29 K, Po ~ P c = 21 
kbar. The curve Tp(P) ~ Tf(P) is almost a straight line, which directly indicates that n = 1 is 
adequate in this case for the description of the P-dependence. The slope for Tps(P) at P = 
is estimated from the condition that its magnitude should not exceed T/o/P c w 0.014 as we have 
assumed that is straight one, so as a result we have —0.014 < p < 0. This ignores the presence 
of a maximum. The available experimental data for ZrZn2 do not give clear indication whether a 
maximum at (T m , P m ) exists. If such a maximum were at P = we would have po = 0, whereas 
a maximum with T m ~ Tps(0) and P m <C Po provides us with an estimated range < po < 0.005. 
The choice po = gives 7 ~ 0.02 and 71 w 0.01, but similar values hold for any |p | < 0.003. 
The multicritical points A and C cannot be distinguished experimentally. Since the experimental 
accuracy [13] is less than ~ 25 mK in the high-P domain (P ~ 20 — 21 kbar), we suppose that 
Tc ~ 10 mK, which corresponds to k ~ 10. We employed these parameters to calculate the T — P 
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diagram using p = and 0.003. The differences obtained in these two cases are negligible, with 
both phase diagrams being in excellent agreement with experiment. 

Phase diagram of ZrZn2 calculated directly from the free energy ([7]) for n = 1, the above mentioned 
values of T s , Po, T/o, and values of 7 fa 0.2 and 71 ~ 0.1 which ensure po ~ is shown in Fig. 
2. Note, that the experimental phase diagram [13] of ZrZn2 looks almost exactly as the diagram 
in Fig. 2, which has been calculated directly from the model (JTj) without any approximations and 
simplifying assumptions. The phase diagram in Fig. 2 has the following coordinates of characteristic 
points: P A ~ P c = 21.42 kbar, P B = 20.79 kbar, P c = 20.98 kbar, T A = T F (P C ) = T FS (P C ) = K, 
T B = 0.0495 K, T c = 0.0259 K, and T F3 (0) = 0.285 K. 

The low-T region is seen in more detail in Fig. 3, where the A, B, C points are shown and the order 
of the FM-FS phase transition changes from second to first order around the critical end-point C. 
The T F s{P) curve, shown by the dotted line in Fig. 3, has a maximum T m = 0.290 K at P = 0.18 
kbar, which is slightly above Tps(0) = 0.285 K. The straight solid line BC in Fig. 3 shows the first 
order FM-FS phase transition which occurs for P B < P < Pc- The solid AC line shows the first 
order N-FS phase transition and the dashed line stands for the N-FM phase transition of second 
order. 

Although the expanded temperature scale in Fig. 3, the difference [T m — T FS (0)] = 5 mK is hard to 
see. To locate the point max exactly at P = one must work with values of 7 and 71 of accuracy up 
to 10 -4 . So, the location of the max for parameters corresponding to ZrZn2 is very sensitive to small 
variations of 7 and 71 around the values 0.2 and 0.1, respectively. Our initial idea was to present a 
diagram with T m = T F s(0) = 0.29 K and po = 0, namely, max exactly located at P = 0, but the final 
phase diagram slightly departs from this picture because of the mentioned sensitivity of the result 
on the values of the interaction parameters 7 and 71. The theoretical phase diagram of ZrZn 2 can be 
deduced in the same way for p — 0.003 and this yields T m = 0.301 K at P m = 6.915 kbar for initial 
values of 7 and 71 which differs from 7 = 271 = 0.2 only by numbers of order 10~ 3 — 10 -4 [18]. This 
result confirms the mentioned sensitivity of the location of the maximum T m towards slight variations 
of the material parameters. Experimental investigations of this low-temperature/low-pressure region 
with higher accuracy may help in locating this maximum with better precision. 

Fig. 4 shows the high-pressure part of the same phase diagram in more details. In this figure the 
first order phase transitions (solid lines BC and AC) are clearly seen. In fact the line AC is quite flat 
but not straight as the line BC. The quite interesting topology of the phase diagram of ZrZn 2 in the 
high-pressure domain (P B < P < Pa) is not seen in the experimental phase diagram [13] because of 
the restricted accuracy of the experiment in this range of temperatures and pressures. 

These results account well for the main features of the experimental behavior [T3], including the 
claimed change in the order of the FM-FS phase transition at relatively high P. Within the present 
model the N-FM transition is of second order up to Pc ~ P c - Moreover, if the experiments are 
reliable in their indication of a first order N-FM transition at much lower P values, the theory can 
accommodate this by a change of sign of 6/, leading to a new tricritical point located at a distinct 
P tr < Pc 011 the N-FM transition line. Since Tc > a direct N-FS phase transition of first order is 
predicted in accord with conclusions from de Haas-van Alphen experiments [H] and some theoretical 
studies [ID]. Such a transition may not occur in other cases where Tc = 0. In SFT (n = 2) the 
diagram topology remains the same but points B and C are slightly shifted to higher P (typically by 
about 0.01 0.001 kbar). 
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8.3 UGe 2 



The experimental data for UGe2 indicate T/o = 52 K, P c = 1.6 GPa (= 16 kbar), T m = 0.75 K, 
P m « 1.15 GPa, and P 0c ps 1.05 GPa [2j |3j SI [5] . Using again the variant n = 1 for T/(P) and the 
above values for T m and P 0c we obtain 7 ps 0.0984 and 71 ~ 0.1678. The temperature T c ~ 0.1 K 
corresponds to k ~ 4. 

Using these initial parameters, together with T s = 0, leads to the T — P diagram of UGE 2 shown 
in Fig. 5. We obtain T A = K, P A = 1.723 GPa, T B = 0.481 K, P B = 1.563 GPa, T c = 0.301 
K, and Pc = 1.591 GPa. Figs. 6 and 7 show the low-temperature and the high-pressure parts of 
this phase diagram, respectively. There is agreement with the main experimental findings, although 
P m corresponding to the maximum (found at ~ 1.44 GPa in Fig. 5) is about 0.3 GPa higher than 
suggested experimentally pfl, [5]. If the experimental plots are accurate in this respect, this difference 
may be attributable to the so-called (T x ) meta- magnetic phase transition in UGe2, which is related 
to an abrupt change of the magnetization in the vicinity of P m . Thus, one may suppose that the 
meta-magnetic effects, which are outside the scope of our current model, significantly affect the shape 
of the T FS (P) curve by lowering P m (along with P B and Pc). It is possible to achieve a lower P m 
value (while leaving T m unchanged), but this has the undesirable effect of modifying P c0 to a value 
that disagrees with experiment. In SFT (n = 2) the multi-critical points are located at slightly 
higher P (by about 0.01 GPa), as for ZrZn 2 . Therefore, the results from the SFT theory are slightly 
worse than the results produced by the usual linear approximation (n = 1) for the parameter t. 

8.4 Two types of ferromagnetic superconductors with spin-triplet elec- 
tron pairing 

The estimates for UGe 2 imply jiK ~ 1.9, so the condition for Tps(P) to have a maximum found 
from Eq. (fT7|) is satisfied. As we discussed for ZrZn 2 , the location of this maximum can be hard 
to fix accurately in experiments. However, P c q can be more easily distinguished, as in the UGe 2 
case. Then we have a well-established quantum (zero-temperature) phase transition of second order, 
i.e., a quantum critical point at some critical pressure Pq c > 0. As shown in Sec. 6, under special 
conditions the quantum critical points could be two: at the lower critical pressure P 0c < P m and 
the upper critical pressure Pq c < P m . This type of behavior in systems with T s = (as UGe 2 ) 
occurs when the criterion (I23p is satisfied. Such systems (which we label as U-type) are essentially 
different from those such as ZrZn 2 where 71 < 7 and hence Tps(0) > 0. In this latter case (Zr- 
type compounds) a maximum T m > may sometimes occur, as discussed earlier. We note that the 
ratio 7/71 reflects a balance effect between the two if)-M interactions. When the trigger interaction 
(typified by 7) prevails, the Zr-type behavior is found where superconductivity exists at P = 0. 
The same ratio can be expressed as Jo/SqMq, which emphasizes that the ground state value of the 
magnetization at P = is also relevant. Alternatively, one may refer to these two basic types of 
spin-triplet ferromagnetic superconductors as "type I" (for example, for the "Zr-type compounds), 
and "type II" - for the U-type compounds. 

As we see from this classification, the two types of spin-triplet ferromagnetic superconductors have 
quite different phase diagram topologies although some fragments have common features. The same 
classification can include systems with T s 7^ but in this case one should use the more general 
criterion (1291) . 
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8.5 Other compounds 

In URhGe, 7/(0) ~ 9.5 K and Tps(0) = 0.25 K and, therefore, as in ZrZn 2 , here the spin-triplet 
superconductivity appears at ambient pressure deeply in the ferromagnetic phase domain [6j [7J [8] . 
Although some similar structural and magnetic features are found in UGe2 the results in Ref. [8] of 
measurements under high pressure show that, unlike the behavior of ZrZn 2 and UGe2, the ferromag- 
netic phase transition temperature Tp(P) ~ Tf(P) has a slow linear increase up to 140 kbar without 
any experimental indications that the N-FM transition line may change its behavior at higher pres- 
sures and show a negative slope in direction of low temperature up to a quantum critical point Tp = 
at some critical pressure P c . Such a behavior of the generic ferromagnetic phase transition temper- 
ature cannot be explained by our initial assumption for the function Tf(P) which was intended to 
explain phase diagrams where the ferromagnetic order is depressed by the pressure and vanishes at 
T = at some critical pressure P c . The Tps{P) line of URhGe shows a clear monotonic negative 
slope to T = at pressures above 15 kbar and the extrapolation [8J of the experimental curve Tps(P) 
tends a quantum critical point Tps{P' ol ?) = at Pq c ~ 25 — 30 kbar. Within the framework of the 
phenomenological theory flBJ this T — P phase diagram can be explained after a modification on 
the Ty(P)-dependence is made, and by introducing a convenient nontrivial pressure dependence of 
the interaction parameter 7. Such modifications of the present theory are possible and follow from 
important physical requirements related with the behavior of the /-band electrons in URhGe. Unlike 
UGe2, where the pressure increases the hybridization of the 5/ electrons with band states lading to a 
suppression of the spontaneous magnetic moment M, in URhGe this effects is followed by a stronger 
effect of enhancement of the exchange coupling due to the same hybridization, and this effect leads 
to the slow but stable linear increase in the function Tp(P) [8J. These effects should be taken into 
account in the modeling the pressure dependence of the parameters of the theory (J7J) when applied 
to URhGe. 

Another ambient pressure FS phase has been observed in experiments with UCoGe j9]. Here the 
experimentally derived slopes of the functions Tp(P) and Tps(P) at relatively small pressures are 
opposite compared to those for URhGe and, hence, the T — P phase diagram of this compound can 
be treated within the present theoretical scheme without substantial modifications. 

Like in UGe2, the FS phase in Ulr [12] is embedded in the high-pressure/low-temperature part of 
the ferromagnetic phase domain near the critical pressure P c which means that Ulr is certainly a 
U-type compound. In UGe 2 there is one metamagnetic phase transition between two ferromagnetic 
phases (FM1 and FM2), in Ulr there are three ferromagnetic phases and the FS phase is located in 
the low-T/high-P domain of the third of them - the phase FM3. There are two metamagnetic-like 
phase transitions: FM1-FM2 transition which is followed by a drastic decrease of the spontaneous 
magnetization when the the lower-pressure phase FM1 transforms to FM2, and a peak of the ac 
susceptibility but lack of observable jump of the magnetization at the second (higher pressure) 
"metamagnetic" phase transition from FM2 to FM3. Unlike the picture for UGe 2 , in Ulr both 
transitions, FM1-FM2 and FM2-FM3 are far from the maximum T m (P m ) so in this case one can 
hardly speculate that the max is produced by the nearby jump of magnetization. Ulr seems to be a 
U-type spin-triplet ferromagnetic superconductor. 
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9 Final remarks 



Finally, even in its simplified form, this theory has been shown to be capable of accounting for a 
wide variety of experimental behavior. A natural extension to the theory is to add a M e term which 
provides a formalism to investigate possible metamagnetic phase transitions [45] and extend some 
first order phase transition lines. Another modification of this theory, with regard to applications 
to other compounds, is to include a P dependence for some of the other GL parameters. The fluc- 
tuation and quantum correlation effects can be considered by the respective field-theoretical action 
of the system, where the order parameters ip and M are not uniform but rather space and time 
dependent. The vortex (spatially non-uniform) phase due to the spontaneous magnetization M is 
another phenomenon which can be investigated by a generalization of the theory by considering 
nonuniform order parameter fields if) and M (see, e.g., Ref. [28J ) . Note that such theoretical treat- 
ments are quite complex and require a number of approximations. As already noted in this paper the 
magnetic fluctuations stimulate first order phase transitions for both finite and zero phase-transition 
temperatures. 
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